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Abstract Illegal logging is a major cause of worldwide 
deforestation, and demands for scientific methods to iden¬ 
tify the geographic origin of timber are increasing. “Den- 
droprovenancing” is one such method, in which the origin 
of unknown wood is estimated by calculating correlations 
of the ring-width series of the unknown wood with refer¬ 
ence trees of known geographic origins. We applied the 
dendroprovenancing method to carbon isotope network 
data of pinyon pines (Pinus edulis and Pinus monophylla) 
from the southwestern United States to test the efficacy of 
using a carbon isotope time series for provenancing wood. 
First, we calculated correlations (t values) between test 
trees temporarily assumed to be of unknown origin and 
reference trees from 13 surrounding sites. Then, we plotted 
the t values on a map. When provenancing was successful, 
the tested trees showed the strongest correlation with refer¬ 
ence trees from sites close to the actual origins of the test 
trees, and the correlations decreased with the distance 
between the original sites of test and reference trees. This 
conical distribution of t values enabled provenancing of 
wood with precision of 114-304 km. Although isotope mea¬ 
surement is more expensive and laborious than ring-width 
measurement, our tests of provenancing pinyon pines in the 
southwestern United States showed a higher success rate 
with carbon isotopes. 
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Introduction 

Deforestation accounts for up to one-third of total anthro¬ 
pogenic carbon dioxide emissions, and illegal logging is one 

't 

of the major causes of worldwide deforestation. However, 
the lack of scientific methods to estimate the geographic 
origin of timber has made the identification of illegally 
logged or traded timber technically difficult. Economical 
damages through illegal timber trade in the world are cur¬ 
rently assessed at 10 billion USD per year. 2 Illegal logging 
is severe in nontropical regions such as the Russian Far East 
and the Baltic States, 3,4 as well as in tropical regions includ¬ 
ing Indonesia. For example, 40% of all wood products 
imported by China are supplied by countries with a high 
incidence of illegal logging, such as Russia and Indonesia. 4 
To prevent illegal logging, consumer countries need to 
ensure traceability of the wood supply from the source to 
the market, thereby preventing “laundering” of illegally 
logged timber before it enters consumer markets. Geo¬ 
graphic origin of wood is the most important aspect of 
traceability. 

Genetic, geochemical, and climatic information pre¬ 
served in wood all have the potential to serve as tools for 
provenancing. In general, identification of tree species at 
the genus level is possible by observing the anatomical 
structure under a microscope. DNA analysis 5 allows identi¬ 
fication of wood at the species level. The origin of timber 
can then be constrained to the distribution of the genus or 
species. Because DNA has natural geographic variation, 
one can also estimate the geographic origin of unknown 
wood by comparing the chloroplast DNA with chloroplast 
DNA from reference trees from different regions. 6 However, 
DNA fingerprinting cannot be applied to plantation trees 
because these DNA fingerprints do not match the natural 
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distribution of DNA. An alternative method to identify the 
geographic origin of wood is thus required. 

Geochemical factors, such as isotope ratios and con¬ 
centrations of inorganic elements, have also been used 
to provenance agricultural products. 7 Isotope ratios of 
strontium in wood 8,9 are known to change geographically. 
However, methods based on stable isotope ratios of 
carbon, nitrogen, and oxygen have proved more effective 
than using inorganic elements for provenancing tropical 
woods. 10 

The ring widths of wood reflect the climatic conditions 
in which the source tree grew, and this knowledge has been 
used for paleoclimate reconstruction. ’ Local climatic fluc¬ 
tuations are reflected in ring-width fluctuations, resulting in 
region-specific ring-width patterns. Trees growing close to 
each other at distances less than a few hundred kilometers 
show higher similarity than the trees growing farther apart. 13 
Ring-width series can thus be used for provenancing 
wood. 14,15 One can estimate the likely origin of unknown 
wood by calculating correlations of the ring-width series of 
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the unknown wood with reference trees of known origins. 
By plotting the correlation of a series from an unknown site 
with each reference chronology on a map, the estimated 
origin can be identified as the location of maximum correla¬ 
tion. 17 This origin identification method, called “dendro¬ 
provenancing,” is mainly used for provenancing historical 
wooden artifacts. 

Stable isotope ratios of tree rings have been used for 
climate reconstruction in a way similar to use of ring-width 
series. 18,19 Isotope ratios reflect the climatic conditions when 
photosynthesis and ring formation took place. 2021 There¬ 
fore, all timber produced anywhere in the world, including 
tropical timber, 2223 has a unique isotopic fingerprint that 
reflects the climatic history of its origin. 18,19 This uniqueness 
suggests the potential use of stable isotope fingerprints for 
provenancing wood. 

In fact, isotope fingerprinting is already in use for esti- 
mating the origins of various plant-based products. “ Both 
oxygen and hydrogen isotope ratios of plants reflect source 
water and leaf transpiration. 29 Isotope ratios of source water 
and leaf transpiration reflect temperature and vapor pres¬ 
sure deficit signals, respectively. Furthermore, carbon 
isotope ratios in C 3 plants reflect the balance between sto- 
matal conductance and photosynthetic rate/ 0 which is dom¬ 
inated by factors such as relative humidity, soil water status, 
summer irradiance, and temperature. In general, using 
oxygen and hydrogen isotopes to provenance agricultural 
products is more successful than using carbon isotopes. 
Despite recent developments in food and forensic fields, to 
date very few studies have so far applied these fingerprint¬ 
ing methods to wood provenancing. If the timber comes 
from nontropical area and has distinct tree rings, then the 
isotopic fingerprint can be easily expressed as a time series 
of yearly stable isotope ratios. 

Application of the “dendroprovenancing” method to 
tree-ring isotope series may further improve the accuracy 
of wood provenancing, because stable isotope ratios of tree 
rings are generally more coherent (i.e., show less within-site 
variation among the individuals) than ring widths. 3133 High- 


frequency fluctuations can be coherent even among sites 
hundreds of kilometers apart. 34 36 

The objective of this study was to evaluate the possibil¬ 
ity of “isotope dendroprovenancing,” i.e., using stable 
isotope series of tree rings to identify the geographic origin 
of wood. As large-scale network data for tree-ring oxygen 
and hydrogen isotopes were not available, we used carbon 
isotope network data from the southwestern United 
States. 37 


Materials and methods 

Carbon isotope network data 

We used tree-ring carbon isotope network data both previ¬ 
ously published 35 37 39 and currently available to the public 
online (http://www.ncdc.noaa.gov/paleo/treering/isotope/). 
The network consists of tree-ring samples from a set of 
14 pinyon pine (Pinus edulis and Pinus monophylla) sites 
[Kane Spring, Alton, Dry Canyon, Lower Colonias, Aztec, 
Cerro Colorado, Ozena, Hawthorne, Mimbres, Owl 
Canyon, north-central Arizona (NC AZ), northeastern 
Arizona (NE AZ), Gate Canyon, and Lamoille] collected 
across six states in the American Southwest (Fig. 1). Site 
elevations ranged from about 1400 to 2500 m, but the sites 
were typically in mountainous terrain with intervening 
topography exceeding 2500 m. In general, the rings of four 
cores from four of the trees were used in the isotopic analy¬ 
sis. After ring-width measurement for cross-dating, rings 
were separated into pentads (e.g., rings formed in 1900- 
1904, 1905-1909, etc.). Pentads from the four cores from 
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Fig. 1 . Locations of the 14 sites and topography in the southwestern 
United States 
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four trees were pooled into composite samples for each 
site except for NC AZ and NE AZ. The pentads of indi¬ 
vidual trees from NC AZ and NE AZ were analyzed 
without pooling. 38 Stable carbon isotope composition (<5 13 C) 
was expressed as follows: 


where A is the carbon isotope discrimination of the ith 
pentad, and A(Holl) is the standardized carbon isotope dis¬ 
crimination, after Hollstein’s transformation. 

In the second method, high-frequency variation is 

O -\ O /IQ 

extracted by calculating first differences ’ ’ as follows: 
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Dendroprovenancing 


(!) A (FD) = A t - A-i (4) 

Finally, in the third method, the difference between the 
ith pentad and the average of the previous and the following 
pentad is calculated as 


We desired to assess the efficacy of isotope dendroprove¬ 
nancing, and one necessary condition for successful den¬ 
droprovenancing is that the actual origin of the unknown 
tree must be surrounded by sites from which reference 
trees were collected. We therefore chose the 6 pooled <5 13 C 
chronologies in the interior of the 14-site network - 1 
pooled chronology each from Kane Spring, Alton, Dry 
Canyon, Aztec, Cerro Colorado, and Gate Canyon 
(Fig. 1) to test as though they were from unknown sites 
(we will hereafter refer to these as test chronologies). 
For the reference network, we used the remaining 6 
marginal chronologies from Fower Colonias, Ozena, 
Hawthorne, Mimbres, Owl Canyon, and Famoille as 
well as the other 7 chronologies from interior sites for a 
total of 13 chronologies. (We will refer to these 13 chro¬ 
nologies as reference chronologies. Note that there are a 
total of 14 chronologies. For each test chronology, the 
remaining 13 chronologies compose a specific reference 
network.) 

By calculating correlations of a test chronology with each 
of the other 13 reference chronologies, the geographic 
origin of the test chronology can be estimated by the den¬ 
droprovenancing method. 13,14,16,17 Both ring widths and <5 13 C 
capture high-frequency (less than 10 years) common forcing 
(macroclimate) better than low-frequency (on the order of 
25 years or more) forcing. ’ ’ Therefore, conventional 
dendroprovenancing techniques preferentially utilize high- 
frequency variations, and we reason that this should apply 
to £ 13 C dendroprovenancing as well. 

As a first step in isotope dendroprovenancing, we 
removed the influence of the Suess effect [effects of chang¬ 
ing <5 13 C of atmospheric C0 2 (<5 L, C a )] from tree-ring S l C 
((5 13 C xr ). For this, we used <5 13 C a data to calculate carbon 
isotope discrimination (zi ). 19,40 


<5 13 C a - d) 13 C TR 
1 + <S 13 C T r/1000 



Then, we extracted high-frequency variations from the 
carbon isotope discrimination. We investigated several 
methods to accomplish this. In the first method, the loga¬ 
rithm of the ratio of ring parameter of a given year to the 
following year (Eq. 3), 41 is used for extracting high-fre¬ 
quency variation as follows: 


— 1 
A+l 7 



A i (AV) = A i 


(A-i + A+i) 


2 



We then calculated the correlation coefficient (r) of high- 
frequency variations between the test chronology and each 
reference chronology. We calculated t values for each refer¬ 
ence chronology as follows: 


t = 




n-2 
1 -r 2 



where n is the number of pentads of the test tree that are 
common with pentads of the reference trees. 

To utilize high-frequency variations in ring-width data 
[RWi(AV)], we processed ring-width data in exactly the 
same way as described above. For both these methods, we 
used data from the corresponding tree rings formed during 
1790-1984 (n = 37) for calculating t values of site chronolo¬ 
gies from the six interior sites. This step enabled the direct 
comparison of the dendroprovenancing results obtained 
from these two different parameters. We estimated the 
origin of the test chronology to be the location of the refer¬ 
ence chronology with which the test chronology showed the 
highest t value (as in Figs. 3-5 in Wazny 17 ). 

To compare the signal-to-noise ratios (SNR) of tree-ring 
<5 13 C and width series, 44 we used a two-way analysis of vari¬ 
ance (ANOVA) applied to the high-frequency variations of 
the five site chronologies [A(AV), RW^AV)] within 300 km 
of Kane Spring (Kane Spring, Alton, Dry Canyon, Aztec, 
and Gate Canyon). Again, we used data from the overlap¬ 
ping tree rings formed during 1790-1985 (n = 37). We 
expected correlations (t values) between two site chronolo¬ 
gies to decrease with the distance between the two sites in 

1 Q 

a similar way as Noda observed with ring width. To verify 
this, we calculated and plotted the t values of all possible 
combinations of two chronologies from the 14 sites ( 14 C 2 = 
91) against the distance. 

Twelve of the site chronologies were constructed from 
four pooled trees each, as already described. However, only 
one tree is expected to be available in most wood prove- 
nancing cases. To examine whether using one tree instead 
of a group of four trees as the test tree can give successful 
provenancing results, we also used the carbon isotope time 
series from individual trees to estimate the origin of each 
tree from NC AZ and NE AZ in the same manner as we 
provenanced the site chronologies. For calculation of this 
set of t values, we used data from the tree ring pentads 
formed during 1845-1984 (n = 26). 


A (Holl) = log 
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Year 

Fig. 2. Change of high-frequency isotope variation with distance. Kane 
Spring {solid circles) chronology showed more similarity with the chro¬ 
nology from nearby Dry Canyon {open circles) than with the chronol¬ 
ogy from distant Cerro Colorado {triangles) 


Results 

Dendroprovenancing of carbon isotope and ring-width 
site chronologies 

Provenancing methods based on Eqs. 3-5 all gave similar 
results. However, for both isotope and ring-width series, 
Eq.5 gave the highest t values, and we therefore adopted 
this equation for all results and figures reported in this 
article unless otherwise specified. Carbon isotope chronolo¬ 
gies from nearby regions showed very similar high-fre¬ 
quency variation [4(AV)]. For example, in Fig. 2 we plotted 
4(AV) from Kane Spring and Dry Canyon, which are 
located 150 km apart (see Fig. 1). However, the similarity 
decreased at Cerro Colorado, which is located 347 km away 
from Kane Spring. 

We considered provenancing to be successful when the 
P value of a correlation maximum was less than 2 x 10“ 7 . 
Threshold t values to meet this level of confidence were 6.45 
in Figs. 3 and 4 {n - 37) and 7.18 in Fig. 5 {n - 26). Such 
high t values (>6.45) were observed within 450 km of the 
actual origin of a test chronology, and t values decreased 
with distance from the origin (Fig. 4), resulting in a conical 
rvalue distribution (Fig. 3a,c,e). For example, the £ 13 C (test) 
chronology from Kane Spring showed the highest t value 
with the reference chronology from Dry Canyon, which is 
150 km away. The t values decreased gradually on the 
western side over the Colorado River drainage basin and 
steeply on the eastern side over the Rocky Mountain range 
(Fig. 3a). The origins of test chronologies were successfully 
estimated as the locations of t value maxima (Figs. 3a,c,e, 
5). This finding strongly indicates the efficacy of carbon 
isotope dendroprovenancing. Four of the six pooled c> Ll C 


chronologies tested (from Kane Spring, Alton, Dry Canyon, 
and Gate Canyon) gave successful provenancing results; 
and the remaining two chronologies (from Aztec and Cerro 
Colorado) gave unsuccessful provenancing results; these 
were “complacent” (low interpentadal variation) chronolo¬ 
gies, and they showed low t values of less than 5.2 with all 
the reference chronologies (Figs. 3g, 4). Among the combi¬ 
nations of two individual carbon isotope chronologies that 
showed successful provenancing {t > 6.45), averaged t values 
were 7.29, 7.26, and 7.74 for Eqs. 3, 4, and 5, respectively. 

When we used ring width instead of <5 13 C for dendroprov¬ 
enancing, the results were less successful. For all test chro¬ 
nologies except for the chronology from Aztec, we failed to 
observe the conical t value distributions with the maxima 
greater than 6.45 (Figs. 3b,d,f, 4). For the test chronology 
from Aztec, the highest t value {t - 7.2; Fig. 4) was observed 
at NC AZ, which is located 522 km away from Aztec (Fig. 
3h). There were eight sites closer to Aztec than NC AZ; 
however, they all showed lower t values than NC AZ. Thus, 
for our study in the southwestern United States, carbon 
isotope dendroprovenancing was more successful than ring- 
width dendroprovenancing. 

Decrease of t values with the distance 

The t values tended to decrease with distance from the 
origin of the test chronologies (see Fig. 4). This trend was 
more evident with carbon isotope ratios than with ring 
widths. We observed high t values exceeding 7 with carbon 
isotopes at distances of 150-300 km, showing the efficacy of 
using tree-ring c) 13 C for provenancing wood. However, 10 of 
14 such combinations of chronologies within 300 km of each 
other showed t values lower than 5; this tended to happen 
especially when high mountains lay between the two chro¬ 
nologies. With ring widths, we observed t value maxima 
exceeding 6.45 at distances 198-522 km away from the 
actual origin of test chronologies (Fig. 4). In other words, 
maxima for ring-width chronologies were dispersed over 
wider areas (-522 km) than maxima for carbon isotope 
chronologies (-304 km). Widely dispersed t value maxima 
for ring width are also reported with European oaks, ’ , 
which is problematic because dispersed distribution hinders 
precise determination of geographic origin. 

Signal-to-noise ratios (SNR) of high-frequency varia¬ 
tions of carbon isotope discrimination [Zl 13 i (AV)] and ring 
width [RWi(AV)] were 8.1 and 4.6, respectively. The higher 
SNR of 4 3 i(AV) that we observed in our study confirms 
earlier findings that tree-ring (5 Ll C is more coherent than 
ring width. 3133 

Isotope dendroprovenancing of individual trees 

At first, we used 14 pooled + 8 individual = 22 series of 
4(AV) for correlation calculation. Four of eight individual 
trees from NC AZ and NE AZ showed t value maxima 
greater than 7.18. However, the other four trees did not 
show t values greater than 7.18 with any of the reference 
series. Because high t values exceeding 7.18 (P < 2 x 10” 7 ) 
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Fig. 3. Results of isotope and 
ring-width dendroprovenancing 
of test chronologies from Kane 
Spring (a, b), Alton (c, d), Dry 
Canyon (e, f), and Aztec (g, h). 
The size of the circles of the 
reference source sites in each 
panel represents the strength of 
the relationship (t value) with 
the site used as the unknown 
source. The color contour plot 
in each panel illustrates t value 
distribution 
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are required for successful provenancing, the 4 low-sensitiv¬ 
ity individual series from NC AZ and NE AZ were excluded 
from the analysis, leaving 18 reference series (=14 pooled + 
4 individual). 

As expected, sensitivity of high-frequency variation 
[^(AV)] to environmental forcing was different among the 
eight individuals. Nevertheless, as few as 28 c) lj C analysis 
data points per individual tree were sufficient to pinpoint 
the geographic origin of the four high-sensitivity series to 
an accuracy of 114-258 km (see Fig. 5). The rate of success 
for provenancing individual carbon isotope series was 50% 
(4 of 8) in our case. The rate of success for provenancing 
individual ring-width series of European oak is reported to 
be 8.9%. 16 


Discussion 

Comparison between isotope and ring-width 
dendroprovenancing 

Tree-ring parameters, such as ring widths or stable isotope 
ratios, reflect common environmental forcing as well as 
random factors within a site. Environmental forcing and 
random factors can be defined as signal and noise, respec¬ 
tively, and we believe our observations of higher SNR of 
4(AV) compared to RWj(AV) indicate less influence of 
random factors on carbon isotope discrimination than on 
ring width. 
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Fig. 4. Decrease of correlations 
(t values) with the distance 
between two chronologies. 

There are 14 C 2 = 91 possible 
combinations of 2 sites of 14 
sites, and all the t values 
generated are shown. The 
origins of test trees could be 
estimated at 150-304 km 
accuracy (, solid-line box) using 
carbon isotopes and 198-522 km 
accuracy ( dashed-line box) with 
ring widths 
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Fig. 5. Results of isotope 
dendroprovenancing of four 
single trees from Arizona 
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Dendroclimatology with ring width is successful when 
the trees growing in stressful environments, such as high 
altitudes or arid environments, are used. 11,12 However, for¬ 
estry is normally practiced at low altitudes where the envi¬ 
ronment is optimal for tree growth. Most of the wood we 
use, whether modern or historical, comes from low alti¬ 
tudes. In such environments, the difference in coherence 
between tree-ring <5 13 C and ring width is expected to become 
wider. For example, isotope chronologies of the trees from 
low altitudes in a humid climate, such as oaks from east 
England and redwood and pine from coastal California and 
Oregon, are known to show better coherence than ring- 
width chronologies. 31 " 33 Even in semiarid Arizona, where 
the difference between isotope ratio and ring width is 
expected to become smaller, we observed a greater success 
rate for dendroprovenancing based on <5 l3 C versus ring 
width. This precision is expected to offer a great advantage 
for provenancing of timber or wooden artifacts. Further¬ 


more, isotope dendroprovenancing requires fewer trees per 
site for building reference networks, 45 and single trees can 
be provenanced with a higher success rate than methods 
based on ring-width chronologies. Disadvantages are that 
isotope measurement is more time consuming (~8 min for 
one measurement) and costly (approximately $5-10 US 
worth of consumables spent for one measurement) than 
ring width measurement (less than 10 s for one measure¬ 
ment at almost no cost). 

Influence of distance and topography between the two 
sites on t values 

As clearly seen in our results, t values are not a simple func¬ 
tion of the distance between the two given sites. We saw a 
strong influence of intersite topography (Fig. 1) on t values 
as evidenced by the gradual decrease in t value moving 
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away from Kane Spring over the Colorado River drainage 
basin compared to the steep decrease over the Rocky 
Mountain range (Fig. 3a). Among the 91 possible combina¬ 
tions of 2 of the total 14 sites, there were 14 combinations 
in which the distance between the 2 sites was less than 
300 km (Fig. 4). Four of these 14 combinations showed suc¬ 
cessful provenancing results (t > 6.45: Alton-Kane Spring, 
Dry Canyon-Kane Spring, Gate Canyon-Kane Spring, NC 
AZ-NE AZ), but the other 10 did not show t values signifi¬ 
cantly high enough to enable provenancing. 

In eight of the ten unsuccessful combinations, high 
mountains exceeding 2500 m lie between the sites, poten¬ 
tially lowering the t values (Aztec-Dry Canyon, Aztec- 
Lower Colonias, Cerro Colorado-Mimbres, Cerro 
Colorado-NE AZ, NE AZ-Mimbres, Alton-NC AZ, Cerro 
Colorado-Dry Canyon, Cerro Colorado-Lower Colonias). 
The rugged topography of Colorado causes large variations 
in climate within short distances, and we must keep in mind 
that few climatic generalizations apply to the whole area. 46 

Another reason for unsuccessful provenancing (Aztec- 
Cerro Colorado, Aztec-Dry Canyon, Aztec-Kane Spring, 
Aztec-Lower Colonias) could be the difference in climate 
sensitivity among trees. 39 For example, Aztec A { ( AV) chro¬ 
nology did not show significantly high t values (t > 6.45) with 
any of the 13 reference chronologies (Fig. 3g). The four 
low-sensitivity series from NC AZ and NE AZ also did not 
show a significantly high t value with any of the reference 
chronologies. We believe that selecting climate-sensitive 
trees as references and then correcting for intertree differ¬ 
ence in climate sensitivity may further improve the prove¬ 
nancing accuracy. 

Error was smaller for provenancing results based on c) 13 C 
(150-304 km) than results based on ring width (198-522 km; 
Fig. 4). When c) 13 C provenancing was successful, the test 
chronologies had a maximum t value higher than the thresh¬ 
old t value at one of the neighboring sites close to the actual 
origin. As our experiment did not include S J C chronologies 
from sites less than 104 km apart (Fig. 4), we cannot say 
whether t values increase further at distances less than 
104 km. The highest achievable spatial resolution of isotope 
provenancing methods depends on how much t values 
increase as the distance between test and reference sites 
approaches zero. In our study, we found that provenancing 
based on ring width often gave the highest t values at distant 
sites instead of sites neighboring the actual origin of test 
chronologies. We can therefore say that spatial accuracy of 
provenancing based on ring width was lower than prove¬ 
nancing based on (5 Ll C. One potential reason for this differ¬ 
ence is that (5 Ll C and ring widths are affected by partly 
different climatic parameters. 1110 For example, precipitation 
events are more localized phenomena than temperature- 
related events. Therefore, temperature observation data 
from a given station are correlated with those of surround¬ 
ing stations in wider areas than when precipitation data are 
used. Carbon isotope chronology from eastern Siberia is 
also correlated with observed July temperature data from 
wider areas than July precipitation data (fig. 5D,E in Kirdy- 
anov et al. 47 ). Carbon isotope ratio is known as a good 
indicator of plant water-use efficiency, 20,30 and it is strongly 


influenced by precipitation. In fact, carbon isotope chro¬ 
nologies from the southwestern United States were used for 
reconstructing drought history. 37,48 

Future prospects of isotope dendroprovenancing 

The carbon isotope chronologies we used were pentads, and 
therefore we could not use cycles shorter than 5 years. 
Higher coherence of high-frequency versus low-frequency 
variations has been previously reported for both ring-width 

Q-1 O/' OQ 

and carbon isotope chronologies. ’ ’ However, suppressed 
trees can occasionally experience sudden decreases in 
neighbor effects as a consequence of tree falls, resulting in 
individual-specific variations, especially in the low-fre¬ 
quency domain. Such conditions may explain why dendro- 
provenancing works best with high-frequency variations. 
Annual 23,33,45 and subannual 23,46,49 isotopic variations are 
reported to be coherent between trees at the same site and 
even between distant sites. Therefore, we believe methods 
based on both annual and subannual isotope variations 
have potential in dendroprovenancing and may even have 
a higher success rate than the methods based on pentads 
reported in this article. In theory, the time-resolution of 
tree-ring isotope archives can reach the monthly level at 
best. 50 

Both SNR and t value maxima of A { ( AV) were higher 
than those of RW^AV). The higher is the SNR of the tree¬ 
ring parameter (or the higher is the correlation coefficient), 
the higher the t value becomes. Because the signal ampli¬ 
tude (macroclimate) is common among the trees and con¬ 
stant within a site, we believe that reduction of noise is the 
key to improve accuracy of isotope dendroprovenancing. 
The noise of 4(AV) is primarily caused by within-site vari¬ 
ability of tree-ring isotope ratios. 23,33,49,51 For example, SNR 
calculated from annual c> Ll C and <5 18 O series from five 
Quercus crispula individuals from the same site are 7.8 and 
14.6, respectively. The noise is also affected by the analyti¬ 
cal error in isotope measurement. 

We suggest the following measures to further improve 
the accuracy of isotope dendroprovenancing. 

o/" 

1. Usage of tree-ring oxygen isotopes 1 and/or hydrogen 
isotopes in combination with carbon isotopes to form 
network data. 

2. Selection of reference trees with high sensitivity to 
common forcing. Correction for the difference in sensi¬ 
tivity among the reference trees may further improve 
provenancing results. However, this necessitates isotopic 
analysis of individual trees and excludes the option of 
pooling. 

3. Usage of high-frequency isotope variation shorter than 
5 years. Annual or even subannual fluctuations are 
coherent among trees within a site, 23,33,48,49 and therefore 
dendroprovenancing may work with shorter frequencies 
than the pentads used in our study. This point is espe¬ 
cially important when the number of available tree rings 
for provenancing is limited. 

4. Development of reference isotope chronologies with 
higher spatial resolution. If future investigations reveal 
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that t values continue to increase at distances less than 
104 km (see Fig. 4), we may be able to provenance trees 
with higher spatial accuracy. 

We believe our study demonstrates the potential of iden¬ 
tifying timber with falsely labeled origins at a spatial accu¬ 
racy of less than 300 km through development of a tree-ring 
isotope network in areas affected by illegal logging, such as 
the Russian Far East and Baltic States. By analyzing disk 
or core samples with exact information about geographic 
origin (longitude and latitude), an isotope lab with sufficient 
facilities, staff, and funding should be able to develop such 
a network for one illegally logged region within a few years. 
However, for this method to be applied practically to market 
timbers, improvement of isotope analysis techniques, such 
as laser-ablation isotope ratio mass spectrometry, 52 a method 
for rapid tree-ring isotope analysis at subannual resolution, 
is necessary. In the tropics, isotope dendroprovenancing 
may be first applied to tree species that form annual rings, 
such as teak. 22 2 ' Teaks are one of the most illegally logged 
species in Indonesia and Burma. 53 Thus, there is a strong 
need for provenancing teak timber. In the future, laser- 
ablation IRMS may also hold the key for provenancing 
tropical timber that lacks tree rings, because even visually 
ringless tropical timber is reported to have annual isotope 
rings. 23 
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